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Preface 

In  this  thesis,  I have  proposed  a system  for  in-flight  fixed  wing 
recovery  of  remotely  piloted  vehicles  (RPV).  The  concept  of  in-flight 
recovery  of  one  aircraft  by  another  is  not  original.  This  concept  dates 
back  to  the  1930’s  when  the  Curtis  F9C-2  Sparrowhawk  was  operated  from  a 
trapeze  attached  to  the  U.S.S.  Akron  dirigible.  In  the  late  1940’s  and 
early  1950’s,  two  parasite  fighter  projects  resulted  in  a few  dockings 
with  another  aircraft.  In  developing  the  proposed  in-flight  recovery 
system  for  RPV’s,  I have  drawn  upon  these  parasite  fighter  programs. 

To  analyze  the  aerodynamic  interference  between  the  capturing  air- 
craft and  the  RPV,  I developed  a horseshoe  vortex  model  to  represent  the 
capturing  aircraft.  The  development  of  this  model  follows  the  develop- 
, ment  presented  by  Glauert.  With  the  aid  of  the  horseshoe  vortex  model, 

I calculated  the  static  change  in  lift,  drag,  and  pitching  moment  on  the 

4 

RPV  at  various  points  along  a flight  path  to  the  docking  point. 

I thank  Professor  Harold  C.  Larsen,  Director,  Air  Force  Institute 
of  Technology  Aerospace  Design  Center,  for  his  guidance  and  assistance 
during  the  preparation  of  this  thesis.  This  thesis  is  dedicated  to  my 
wife  for  her  understanding  and  support  during  its  preparation. 


i 


ii 


WJ*  '-.WISWW . mi ;. 


i?rf" 


1>1W!H  t|,.U 


GAW/AE/74-2 


4 

0 

Contents 


Preface  

List  of  Figures  

List  of  Symbols  

Abstract 

I.  Introduction  .....  

II.  Design  Concepts.  

Fixed  Wing  RPV  

Oblique  Wing  RPV  

Capture  System  . 

Skyhook 

Trapeze 

Aft  Support,  Fixed  Wing  RPV 

Aft  Support,  Oblique  Wing  RPV 

III.  Horseshoe  Vortex  Model  of  Mother ship  

Contribution  of  Bound  Vortex  

Contribution  from  Trailing  Vortices 

Calculation  of  ^ 

Calculation  of  q-j 

Total  Induced  Velocity  

Analysis  of  Induced  Velocity  

IV.  Steady  State  Analysis  of  Aerodynamic  Interference 

Variation  of  Lift 

Variation  of  Drag 

Variation  of  Pitch  Stiffness  

Analysis  of  Hypothetical  RPV  

V.  Conclusions  and  Recommendations 

Conclusions 

Recommendations.  ........  

Cited  References 

Appendix  A:  Computer  Program  

Vita 


iii 


4 


jijijj]  iji.jjm 


GAW/AE/74-2 


List  of  Figures 

Page 


Figure 

1 Oblique-Winged  RPV  5 

2 RPV  Capture  System  (Bomb  Bay  Doors) 7 

3 RPV  Capture  System  (Aft  Doors) 8 

4 Fixed-Wing  RPV  with  Skyhook  Extended  10 

5 Skyhook 11 

6 Bomb  Bay  Trapeze 13 

7 Aft  Door  Trapeze 14 

8 Mothership  Horseshoe  Vortex 17 

9 Induced  Velocity  Distribution,  Um  ■ 300  Knots 25 

. 10  Induced  Velocity  Distribution,  * 200  Knots 26 

11  RPV  Flight  Path 30 

12  Steady  Level  Flight  of  RPV 32 

13  Definition  of  Angles  . . 33 

14  Definition  of  Pitch  Moment  Parameters 38 

15  Velocity  Change  of  RPV 44 

16  Lift  Change  of  RPV 45 

17  Change  in  RPV  Coefficient  of  Lift 46 

18  Angle  of  Attack  Change  .....  47 

19  Change  in  Coefficient  of  Drag 48 

20  Drag  Change  of  RPV 49 

21  Thrust  Change  of  RPV  , 50 

22  Change  in  RPV  Static  Pitch  Stability 52 

23  Closing  Velocity  53 

24  Time  to  Dock 54 

25  Simplified  Flow  Diagram 59 


1 


i 


iv 


H 


GAW/AE/74-2 


List  of  Symbols 


Lift  Curve  Slope 
Aspect  Ratio 
Mean  Aerodynamic  Chord 
Center  of  Gravity 
Coefficient  of  Drag 
Coefficient  of  Drag  at  Zero  Lift 


Coefficient  of  Lift 


Moment  Coefficient 


Derivative  of  Static  Pitch  Moment  with  Respect  to  a 


Oswald  Airplane  Efficiency  Factor 
Location  of  c.g.  (Fig.  14,  p.  38) 

Location  of  wing  mean  aerodynamic  center  (Fig.  14,  p.  38) 

Tail  Incidence  Angle 

Correction  Factor  in  Coefficient  of  Drag 

Distance  Between  c.g.  and  Tail  Aerodynamic  Center 

Distance  Between  Wing  Aerodynamic  Center  and  Tail  Aero- 
dynamic Center 


P(X1,Y1,Z1) 


rl,r2*r3 


Moment 


Point  of  Analysis  (Location  of  RPV) 
Induced  Velocity  from  Bound  Vortex 
Induced  Velocity  from  +s'  trailing  Vortex 
Induced  Velocity  from  -s'  Trailing  Vortex 
Vectors  Defined  in  Fig.  8,  p.  17 
Ratio  s'/s 


v 


I 

k 


I 


[ 

i 


l 

* 

\ 

f 


* 


GAW/AE/74-2 


8 

S' 

S 

z 

zt 


Subscripts 

ac 

t 

W 


* Greek  Letters 
a 

aLo 

a 

00 

C 


T 


* 

0 

Semi-span  of  Wing 

Semi-span  of  Trailing  Vortex  Filaments 
Wing  Area 

Wing  Offset  from  c.g.  (Fig.  14,  p.  38) 
Tail  offset  from  c.g.  (Fig.  14,  p.  38) 


Aerodynamic  Center 

Tail 

Wing 

Free  Stream  Condition 


Angle  of  Attack 

Angle  of  Attack  Zero  Lift 

Free  Stream  Angle  of  Attack 

Downwash  Angle  from  RPV  Wing 

Circulation  Around  Mid-span  of  Wing 

Elliptic  Variation  of  Circulation  on  Wing  Span 

Density  at  Altitude 

Density  at  Sea  Level 

Angle  Between  Relative  Velocity  and  Free  Stream 
Velocity  (Fig.  13,  p.  33) 


GAW/AE/74-2  f 

4 | 

# r 

Abstract  H 

A potential  approach  for  recovery  of  remotely  piloted  vehicles  (RPV) 
is  the  in-flight  docking  of  the  RPV  with  a mothership.  Some  design 
concepts  are  presented  for  in-flight  recovery  systems.  A steady  aero- 

j 

dynamic  analysis  waj  performed  on  a hypothetical  RPV  as  it  approached 

the  mothership  for  docking.  The  correct  trend  in  the  variation  of  lift  ; 

i 

and  static  pitch  moment  of  the  RPV  was  predicted  using  a horseshoe  vortex 
model  to  represent  the  mothership. 
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I.  Introduction 

Recent  articles  in  the  open  literature  indicate  numerous  advantages 
of  using  Remotely  Piloted  Vehicles  (RPV's)  in  place  of  piloted  vehicles. 
One  of  the  main  advantages  of  using  RPV's  is  that  they  can  operate  in  a 
desired  area  without  potential  loss  of  a pilot.  One  of  the  major  dis- 
advantages in  using  RPV's  is  the  number  of  personnel  and  amount  of  equip- 
ment presently  required  to  launch  and  recover  a single  RPV.  During  the 
present  Department  of  Defense  budget  limitations,  total  funds  available 
are  remaining  relatively  constant  while  the  manpower  costs  are  rising 
rapidly.  By  reducing  the  number  of  personnel  required  to  recover  RPV's, 
the  cost  of  RPV  operations  may  be  decreased. 

For  the  RPV's  that  are  air  launched  from  a "mother"  aircraft,  it 
may  be  more  cost  effective  to  have  the  launch  aircraft  also  act  as  the 
recovery  system.  Under  current  practices  the  launch  and  recovery  systems 
are  completely  separate  and  distinct.  The  development  of  a integrated 
launch/recovery  system  provides  the  potential  for  a cost  savings  by  re- 
ducing the  number  and  type  of  land,  air,  or  water  vehicles  required  to 
carry  out  RPV  operations. 

One  possible  concept  for  an  integrated  launch /recovery  system  is  in- 
flight capture  of  the  RPV  by  a fixed  wing  aircraft.  Mid-air  capture  of 
one  aircraft  by  another  aircraft  is  not  a new  concept.  Both  the  XF-85  and 
the  XF-84  projects  of  the  late  1940' s and  early  1950' s showed  that  it 
was  physically  possible  for  one  aircraft  to  rendezvous  and  dock  with 
another  aircraft.  However,  both  programs  showed  that  mid-air  rendezvous 
and  docking  was  not  simple.  As  the  docking  aircraft  approached  the  cap- 
ture mechanism  of  the  mother  aircraft,  the  aerodynamic  interference  be- 
tween the  two  aircraft  drastically  increased  the  occurance  of  adverse 
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effects  on  the  docking  aircraft  which  was  much  smaller  than  the  mother 
aircraft.  The  drag  on  the  docking  aircraft  increased  making  it  difficult 
to  maintain  a steady  closing  speed.  A degradation  in  the  stability  of 
the  docking  aircraft  occurred  as  the  capture  mechanism  was  approached. 

Before  flight  tests  are  performed  on  potential  launch/recovery 
systems,  an  analysis  of  the  aerodynamic  interference  between  the  aircraft 
involved  should  be  performed.  Capture  wind  tunnel  tests  should  be  per- 
formed to  verify  the  analysis  or  provide  data  for  the  analysis.  The  re- 
design necessary  as  a result  of  flight  testing  may  be  minimized  by  per- 
forming an  analysis  and  wind  tunnel  test  of  the  aerodynamic  interference 
before  the  aircraft  is  built  or  modified. 

The  objective  of  this  thesis  is  to  perform  a steady  state  aerodynam- 
ic analysis  of  the  interference  between  the  two  aircraft  of  a proposed 
integrated  launch/recovery  system  for  RPV's.  In  this  thesis,  seme  de- 
sign concepts  for  the  equipment  necessary  for  in-flight  docking  of  a 
fixed  wing  and  oblique  wing  RPV  are  proposed.  In  the  analysis  of  the 
aerodynamic  interference  the  mothership  wing  is  modeled  by  a horseshoe 
vortex  system  with  elliptical  distribution  of  lift  on  the  wing.  The 
fuselage  and  tail  effects  on  the  mothership  on  the  RPV  are  neglected  in 
the  analysis. 

The  analysis  of  the  interference  between  the  two  aircraft  consists 
of  computing  the  steady  state  values  of  RPV  lift,  drag,  and  aerodynamic 
pitch  moment  at  positions  relative  to  the  mothership  as  the  point  of 
analysis  is  moved  to  the  docking  point.  The  thesis  is  divided  into  three 
major  areas: 

a.  Design  concepts 

b.  Horseshoe  vortex  model  of  mothership 


c.  Steady  state  analysis  of  aerodynamic  interference 
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*11.  Design  Concepts 


Aerial  capture  of  one  aircraft  by  another  is  not  a new  concept.  The 
author  has  found  in  the  literature  three  systems  which  were  designed  for 
and  performed  aerial  capture.  The  first  system,  the  Curtiss  F9C-2 
Sparrowhawk,  was  operated  from  a trapeze  installation  aboard  the  U.S.S. 
Akron  dirigible  in  1933.  The  second  system  was  the  XP--J5  "GOBLIN"  de- 
signed by  McDonnell  Aircraft  Corporation  (1944-1949).  The  XP-85  was  de- 
signed as  a parisite  fighter  to  be  carried  internal  to  the  B-36.  Only 
two  XP-85  prototype  were  built  and  flight  tested.  The  third  system  was 
the  F-84F  FICON,  also  designed  to  be  carried  by  the  B-36.  The  FICON  was 
a modified  Republic  F-84  which  was  flight  tested  in  March  1953.  In 
addition,  F-84  aircraft  were  designed  to  be  wing  tip  carried  on  the  B-50, 
but  successful  "hook-on"  of  two  aircraft  was  never  made. 

In  the  discussion  of  design  concepts  for  in-flight  recovery  of  RPV's, 
two  different  types  of  aircraft  will  be  discussed.  The  first  type  of 
RPV  is  a fixed  wing  aircraft  and  the  second  is  an  oblique  wing  RPV  on 
which  the  wing  pivots  as  a unit. 


Fixed  Wing  RPV 

The  fixed  wing  RPV  recovery  is  based  on  aircraft  like  the  Firebee 
drone.  The  Firebee  drone  aircraft  has  the  ving  fixed  to  the  fuselage. 

For  fixed  wing  RPV's,  some  mechanism  will  most  likely  have  to  be  added 
to  the  wings  to  allow  them  to  fold  for  storage  inside  the  mothershir. 

The  F-85  project  used  a folding  mechanism  on  its  wings  so  that  it  could 
be  stored  in  the  bomb  bay  of  the  B-36. 

Oblique  Wing  RPV 

Discussion  of  oblique  wing  RPV  recovery  is  based  on  an  aircraft  that 
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has  its  wing  mounted  on  a pivot.  The  pivoted  wing  allows  the  complete 
wing  to  be  rotated  so  that  it  is  parallel  to  the  fuselage,  thus  decreasing 
the  wing  span  for  storage  inside  the  mothership. 

Some  of  the  advantages  for  using  oblique  winged  aircraft  are  presented 
in  an  article  by  Robert  T.  Jones  (Ref  1:66-70).  Jones  proposed  a oblique- 
winged aircraft  as  a candidate  for  the  Super  Sonic  Transport.  The  oblique- 
wing of  his  proposed  aircraft  pivots  as  a unit  to  different  oblique  angles 
corresponding  to  flight  at  different  Mach  numbers.  Jones  gives  the  follow- 
ing advantages  for  an  oblique-winged  aircraft: 

a.  The  maximum  lift-drag  ratio  for  various  Mach  number  ranges  are 
"significantly  higher  than  those  previously  obtained  with  arrow  wing  or 

delta  wing  body  combinations."  * 

b.  The  oblique  elliptic  wing  has  the  optimum  supersonic  wing  plan- 
form  indicated  by  the  drag  reversibility  theorem  of  Von  Karman  and  Hayes. 

c.  The  wing  structure  is  continuous  across  the  pivot  producing  only 
a tension  load  on  the  pivot  with  none  of  the  bending  loads  of  the  "swing 
wing." 

d.  Turning  the  wing  as  a whole  does  not  cause  a rearward  shift  of 
the  center  of  lift  because  the  centroid  of  area  does  not  move  relative 
to  the  center  of  gravity. 

e.  The  optimum  fuselage  shape  for  the  oblique-wing,  which  conforms 
to  the  "area  rule,"  is  nearly  cylindrical  as  compared  to  the  localized 
deep  indentation  of  the  fuselage  for  the  swept  wing. 

An  oblique-winged  aircraft  of  the  type  discussed  by  Jones  is  pro- 
posed for  use  with  the  integrated  launch/recovery  system  discussed  in 
this  thesis  (Fig.  1,  p.  5).  The  pivoted  oblique-wing  will  allow  the 
RPV  to  fly  with  the  wing  straight  for  capture  at  subsonic  speed  and  to 
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fly  at  various  Mach  numbers  while  performing  its  mission.  With  the  RPV 
locked  to  the  mothership  by  the  supports,  the  wing  can  be  pivoted  to  de- 
crease the  span  for  storage  inside  the  mothership  fuselage. 

Capture  System 

Potential  in-flight  capture  systems  can  be  based  on  the  following 
approaches:  a.  launch  and  recovery  of  the  RPV  from  the  wing  of  the 

mothership  with  external  carry,  b.  launch  and  recovery  of  the  RPV  from 
the  tail  of  the  mothership  with  external  carry,  and  c.  launch  and  recov- 
ery of  the  RPV  from  the  fuselage  of  the  mothership  with  internal  carry. 

The  latter  approach  is  the  basis  for  the  concepts  presented  in  this  thesis. 

The  capture  system  is  divided  into  the  trapeze  (carried  on  the  mother 
ship  aircraft)  and  the  skyhook  (carried  on  the  RPV).  The  terms  trapeze 
and  skyhook  were  used  in  the  XP-85  project  and  are  carried  over  by  this 
author.  Two  capture  systems  are  considered  in  this  thesis.  The  first 
system  involves  recovery  through  bomb  bay  doors  in  the  fuselage  (Fig.  2, 
p.  7 ).  The  second  system  involves  recovery  through  a tail  door  (Fig.  3, 
p.  8)*  In  both  capture  systems  the  RPV  flies  from  the  rear  toward  the 
mothership.  During  the  capture  sequence,  the  RPV  approaches  the  trapeze 
(docking  point)  with  the  skyhook  extended.  As  the  RPV  closes  on  the 
trapeze  the  closing  rate  is  reduced.  The  skyhook  is  guided  into  the 
capture  bar  by  two  arms  (see  Fig.  6,  p.  13).  Once  the  RPV  locks  on  to 
the  bar  the  engine  thrust  is  reduced  md  the  RPV  is  then  towed.  While 
the  RPV  is  in  tow,  the  aft  support  is  lowered  and  locked  in  place  to 
stabilize  the  RPV  for  raising  into  the  mothership  for  storage.  Once  the 
aft  supports  are  secured  the  engine  is  shut  off. 
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Figure  2.  RPV  Capture  System  (Bomb  Bay  Doors) 


Figure  3.  RPV  Capture  System  (Aft  Doors) 
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The  skyhook  end  Che  associated  trapeze  is  the  same  for  both  capture 
systems  (see  Skyhook,  p.  9 and  Trapeze,  p.  12  ).  The  aft  support  used 
with  the  boiul>  lay  doors  is  shown  in  Fig.  6 , p.  13  . The  wing  support 
used  with  the  tail  door  is  shown  in  Fig.  7,  p.  14  . 

Skyhook.  The  skyhook  for  the  XP-85  was  located  2 to  7 inches  ahead 
of  and  above  the  center  of  gravity  of  the  aircraft.  During  the  flight 
testing  of  the  XP-85,  with  a trapeze  attached  to  a B-29,  McDonnell 
Aircraft  Corporation  found  that  this  position  was  unsuitable  because  the 
XP-85  became  marginally  stable,  direccionally , wnile  attached  to  the 
trapeze.  McDonnell  recommended  that  any  future  designs  should  locate 
the  skyhook  well  ahead  of  the  center  of  gravity  so  that  the  aircraft  is 
towed  after  attachment  to  the  trapeze  (Ref  2:10).  The  XF-84F  which 
followed  the  XF-85  had  the  skyhook  located  close  to  the  nose  of  the  air- 
craft (Ref  3:  A3  ) . Based  on  the  information  gained  from  the  XP-85 
and  XF-84F  projects,  the  skyhook  should  be  located  on  the  rose  of  the 
RPV  as  shown  in  Fig.  4,  p,  10  . The  exact  location  of  the  skyhook  must 
be  determined  by  an  analysis  of  the  RPV  stability  while  attached  to  the 
trapeze. 

The  skyhook  design  concepts  presented  are  similar  to  the  design 
used  in  the  XF-84F  project.  The  skyhook  design  is  common  to  both  capture 
systems.  The  skyhook  (see  Fig.  5,  p.  11  ) consists  of  a hook  rigidly 
attached  to  the  RPV  at  its  base  and  free  to  rotate  in  the  plane  of 
symmetry  of  the  RPV.  The  hook  is  raised  and  lowered  by  a shock  absorbing 
shaft  attached  to  a slide.  The  slide  travels  along  a track  which  is 
fixed  to  the  RPV.  When  the  skyhook  is  completely  extended  the  slide  is 
locked  in  place  for  docking.  The  shock  absorbing  shaft  is  used  to  absorb 
some  of  the  energy  that  must  be  dissipated  during  docking. 
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Figure  4.  Fixed-Wing  RPV  with  Skyhook  Extended 


Figure  5.  Skyhook 
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In  the  hook  is  a latch  which  locks  onto  the  docking  bar  of  the 
trapeze.  At  the  end  of  the  hook  a guide  is  attached  by  shear  pins  to 
increase  the  contact  area  with  the  attachment  bar  of  the  trapeze. 

Trapeze.  During  docking  the  skyhook  locks  onto  the  trapeze  attach- 
ment bar.  The  trapeze  puts  the  RFV  in  tow  by  the  mothership. 

The  trapeze  is  a extendable  shaft,  similar  to  the  KC-135  refueling 
boom,  with  a docking  head  on  the  end  (Fig.  6,  p.  13  and  Fig.  7,  p.  14  ) 
and  is  common  to  both  capture  systems.  The  docking  head  should  be  at 
least  two  fuselage  diameters  of  the  mothership  below  the  fuselage.  This 
length  was  suggested  as  a result  of  the  flight  tests  of  the  XF-85  project 
(Ref  2:  10)  which  indicated  that  the  distance  was  outside  the  turbulent 
boundary  layer.  Located  above  the  docking  head  are  two  moveable  guide 
wings  for  controlling  the  trapeze. 

The  docking  head  is  attached  to  the  tow  boom  by  a pivot  which  will 
allow  the  docking  head  to  swing  in  the  plane  of  symmetry  of  the  mother- 
ship. A shock  absorber  is  attached  at  one  end  to  the  tow  boom  and  at  the 
other  end  to  the  docking  head.  The  shock  absorber  is  used  to  hold  the 
docking  head  vertical  and  to  absorb  some  of  the  energy  of  the  RFV  hook 
up  with  the  docking  head. 

The  docking  head  is  a fork  with  the  two  prongs  pointed  vertically 
downward,  away  from  the  mothership.  Joining  the  tips  of  the  prongs  is 
the  attachment  bar.  Guides  are  attached  to  the  tips  of  the  prongs  to 
increase  the  contact  area  and  guide  the  skyhook  to  the  attachment  bar. 

With  the  skyhook  locked  on  to  the  attachment  bar  of  the  trapeze 
the  thrust  is  reduced.  While  the  RPV  is  being  towed,  the  aft  support  is 
lowered.  Each  of  the  aft  supports  is  discussed  separately. 

Aft  Support,  Fixed  Wing  RPV.  The  aft  support  is  designed  to  provide 


12 


Bomb  Bay  Trapeze 


14 


pitch  and  directional  stability  during  the  storage  sequence.  The  aft 
support  contacts  the  fixed  wing  RPV  on  the  fuselage  behind  the  center 
of  gravity.  Fig.  6,  p.  13  shows  a schematic  drawing  of  the  aft  support. 
The  aft  support  is  located  at  the  end  of  a extendable  boom  like  the  docking 
head.  A set  of  claw  like  arms  are  used  to  encircle  the  aft  part  of  the 
fuselage  and  lock  in  place. 

Aft  Support . Oblique  Wing  RPV.  For  the  oblique-wing  RPV  proposed 
in  this  thesis  a different  aft  support  may  be  used.  The  aft  support 
shown  schematically  in  Fig.  7,  p.  14  , is  a wing  support.  The  extend- 
able rig  moves  on  a rail  attached  to  the  mothership.  When  the  RPV  is  in 
tow,  the  wing  support  table  is  lowered  in  front  of  the  RPV  wing  and  then 
moved  back  to  enclose  the  wing  on  both  sides  of  the  pivot.  The  wing 
clamp  is  then  locked  in  place  for  raising  into  the  fuselage.  For  storage 
the  wing  support  table  and  wing  are  rotated  so  that  the  wing  is  parallel 
to  the  fuselage. 
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III.  Horseshoe  Vortex  Model  of  Mothership 

In  this  thesis  the  mothership  is  modeled  by  an  elliptical  distribu- 
tion of  horseshoe  vortex  elements  on  the  mothership  wing.  The  induced 
velocity  on  the  RPV  during  the  capture  sequence  was  obtained  from  the 
horseshoe  vortex  model.  Because  the  RPV  was  assumed  to  be  much  smaller 
than  the  mothership,  the  effect  of  the  RPV  on  the  mothership  velocity 
environment  was  assumed  to  be  negligible.  Figure  8,  p.17  , shows  the 
coordinate  system  and  the  relationships  used  in  the  horseshoe  vortex 
model  of  the  mothership.  The  induced  velocity  at  point  (XI,  Yl,  Zl) , 
where  point  P represents  the  location  of  the  RPV  wing  from  the  origin  in 
the  mothership,  is  obtained  using  the  Biot-Savart  Law.  The  total  induced 
velocity  q at  the  point  P can  be  obtained  by  summing  the  individual  con- 
tributions of  the  bound  vortex  and  each  of  the  two  trailing  vortex  filaments. 

As  an  Initial  step  in  formulating  the  mothership  model,  an  elliptical 
spanwise  distribution  of  lift  of  the  mothership  wing  was  assumed.  Glauert 
justifies  this  assumption  by  the  following  statements  taken  from  his  book 
(Ref  4:143). 

"The  elliptic  distribution  of  circulation  or  lift  across 
the  span  of  an  aerofoil  is  important,  firstly  because  it  leads 
to  the  minimum  possible  induced  drag  for  a given  total  lilt,  and 
secondly  because  the  load  grading  curves  of  most  aerofoils  of 
conventional  shape  do  not  differ  greatly  from  the  elliptic  form. 

The  results  deduced  from  the  hypothesis  of  elliptic  loading  are 
therefore  the  best  which  can  possibly  occur  and  are  also  a good 
first  approximation  to  those  actually  obtained". 

The  mothership  was  modeled  by  the  wing  only,  no  fuselage  or  tall 
effects  were  considered  in  the  model.  The  wing  was  assumed  to  be  straight 
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Side  View 


Vector  Relationships 
dx  . — d x i 


|P(X1,Y1,Z1) 


NOTE: 


(p  - » XI  i + (Y1  - y)  j + Z1  k 

(p  - OA)  - (XI  - x ) i + (Y1  - sR)  j + Z1  k 

(p  - o&)  - (xi  - x)  i + (y 1 + sr)  j + zi  k 

The  sign  of  the  symbols  used  in  the  vector  expressions 
are  positive  except  for  sR.  Thus  the  value  of  the  symbol  to 

be  substituted  into  the  above  expression  carries  a si<?n. 


Figure  8.  Mothership  Horseshoe  Vortex 
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with  an  elliptic  distribution  of  lift  across  the  span.  The  elliptic 
distribution  of  lift  (elliptic  loading)  is  proportional  to  the  circula- 
tion around  the  airfoil.  The  wing  with  a steady  flow  of  air  passing  over 
it  can  be  represented  by  individual  vortex  lines  which  combine  linearly 
to  produce  the  elliptic  distribution  of  circulation  around  the  airfoil. 
These  individual  vortex  lines  consist  of  bound  vortex  segments,  represent- 
ing the  wing  itself,  and  free  vortex  lines  that  trail  from  the  wings 
trailing  edge.  The  trailing  vortex  lines  form  a trailing  vortex  sheet 
behind  the  wing.  When  the  starting  vortex  is  infinitely  down  stream,  the 
vortex  lines  form  what  is  called  a horseshoe  vortex. 

The  sheet  of  trailing  vortex  lines  is  unstable,  so  downstream  of  the 
airfoil  the  vortex  sheet  rolls  up  into  two  trailing  vortex  lines  which 
have  strength  equal  to  the  sum  of  the  individual  vortex  filaments  that 
make  up  the  trailing  vortex  sheet.  For  an  elliptically  loaded  wing  the 
sum  of  the  individual  horseshoe  vortex  filaments  is  a maximum  at  mid-span 
with  a value  of  T0  and  decreases  elliptically  to  zero  at  the  wing  tips. 

The  trailing  vortices  then  have  a constant  strength  of  rQ.  The  centers 
of  the  trailing  vortices,  for  an  elliptically  loaded  wing,  are  shown  l y 
Milne-Thomson  (Ref  5:206)  to  be  separated  by  a distance  equal  to  it /A 
times  the  wing  span. 

Treating  the  wing  as  a horseshoe  vortex  to  study  the  effects  away 
from  the  wing  itself  is  stated  by  Houghton  and  Brock  (Ref  6:336)  to  give 
good  results  at  distances  greater  than  two  cnord  lengths  from  the  center 
of  pressure  of  the  wing.  The  effect  of  the  mothership,  modeled  by  the 
horseshoe  vortex,  is  found  by  first  determining  the  induced  velocity 
contribution  from  the  bound  vortex,  and  then  finding  the  induced  velocity 
contribution  from  each  of  the  trailing  vortices. 
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Contribution  of  Bound  Vortex 

The  bound  vortex  lies  along  the  y-axis  with  a length  2s,  where  s is 
the  semi-span  of  the  mothership  wing.  The  bound  vortex  induced  velocity 
at  point  P is  called  q^.  For  the  elliptically  loaded  airfoil  the  dis- 
tribution of  circulation  across  the  span  is 

r (y)  - — /s2  - y2  (1) 

s 

where  rQ  is  the  maximum  circulation  at  the  midsection  of  the  wing  (Ref 
7:548).  The  value  of  T0  is  given  by 


2L  2W 

m - - — 

irpV^  S lrpV^  S 


(2) 


(Ref  7:550). 

4 

The  induced  velocity  qj  at  point  P is  then  given  by  the  Biot-Savart 


Law  us 


\ 1 fS  r(y)(dyXri) 

V(r.)  ■ — j — — (cross  product) 

An  -s  |r1|-> 


(3) 


(Ref  7:528). 

Substituting  in  the  expression  for  F(y)  the  induced  velocity  q ^ is 
+ T0  s h2  - y^(dyXr+1) 


-i-  r r 

4ns  -s  lrll 
As  shown  in  Fig.  b f p,  18 f 


(4) 


r1  - (P  - y)  - XI  i + (Y1  - y) j + Z1  k 


(5) 


and  dy  ■ dy  j.  Evaluating  the  cross  product  dy  X r^  gives 

dy  X r^  ■ Z1  dy  L - XI  dy  k.  ((,) 

Also,  |r^|  * /Xl^  + (Y1  - y)^  + Zl^.  Substituting  these  expressions  into 
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the  equation  for  q^  gives' 

a _ i * * 

q -£o_  /*  <9  - y >2  <Z1  dy  1 - XI  dy  k) 

1 4us  -s  (XI2  + (Y1  - y)2  + Zl2)3/2 

The  expression  is  now  factored  into  the  two  components 

; r0Zl  » (s2  - y2)7  dy 

q,  i - J 

x 4ns  -s  (XI2  + (Y1  - y)2  + Zl2)3^2 


and 


«1, 


roX1  js 


(3  - y )?  dy 


4ns  -6  (XI2  + (Y1  - y)2  + ZI2)3/2 


(7) 


(8) 


(9) 


Contribution  from  Trailing  Vortices 

The  trailing  vortex  sheet  formed  by  the  distribution  of  individual 
horseshoe  vortices  across  the  wing  is  unstable  at  the  wing  tips  and  rolls 
up  into  two  trailing  vortex  filament''-.  The  single  trailing  vortex  fila- 
ments that  result  from  the  vortex  sheet  rolling  up,  start  at  points  set 
back  from  t \e  wing  tips.  The  distance  the  trailing  vortex  is  off  set 
from  the  wing  tip  is  a distance  (s  - s')  where  s is  the  wing  semi-span 
and  s'  is  the  semi-span  of  the  trailing  vortex  filaments.  The  two  semi- 
spans s and  s'  can  be  related  by  the  ratio  R * s'/s.  For  an  elliptically 
loaded  wing  Glauert  gives  a ratio  of  s'/s  ■ j (Ref  4:168).  The 
strength  of  each  trailing  vortex  filament  is  the  sum  of  each  element  in 
the  sheet  and  is  equal  to  the  circulation  around  the  center  of  the  wing, 
that  is  rQ. 

Using  the  Biot-Savart Law  with  the  circulation  constant  at  ro  the 
induced  velocity  from  the  +s'  trailing  vortex  and  the  -s'  trailing  vortex 
are  q2  and  q^,  respectively.  To  calculate  ^ and  the  notation  of 
Fig.  8,  p.  17  is  used. 
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Calculation  of  q2»  Beginning  with  the  Biot-Savart  Law  for  a serai- 
infinite  vortex  filament  starting  at  +s',  the  induced  velocity  q2  at  point 
P is 

■+  -*• 
r -•  v ». 

(10) 


r0  (dxXr, 

■ q2(p) " 7 > ~FT^~ 

4ir  o |r2|J 


where  dx*  + dx  i (dx  points  in  the  direction  of  integration  and  is  a 
negative  increment  in  the  integration.) 


(ID 


r2  - (XI  - x)  i + (Y1  - sR)  j + Z1  k , 

and  |r2|  - /(XI  - x)2  + (Y1  - sR)2  + Zl^ 

Evaluation  of  the  cross-product  gives  dx  X r2  * -Zl  dx  j + (Y1  - sR) 

* •+ 
dx  k.  Substituting  these  into  the  integral  equation  for  q2 (p)  the  follow- 

. ing  is  obtained. 

la  21  i - - sl!>^  * (12) 

4l  { [(XI  - x>2  + (XI  - sR)2  + Zl2]3/2 

This  equation  is  broken  into  the  components  and  integrated. 

To  integrate  the  equation  for  q£ (P)  the  denominator  is  expanded  to 

x2  - 2 XI  x + XI2  + (Y1  - sR)2  + Zl2  (13) 

Letting 


a - XI2  + (Y1  - sR)2  + Zl2 


b - -2  XI 
c ■ 1 

and  using  integral  formula  No.  184  from  Ref  8:297; 

t dx  * 2(2cx  + b) 


(14) 


(15) 


where  X ■ a + bx  + cx^ 

2 

q ■ 4ac  - b 


the  following  expression  is  obtained  for  q2^(P)  J' 
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is  opposite  to  that  for  q£.  The  induced  velocity  is  thus 

„»  m _ -fo  f"  <d*  x »3> 

,d<P)  17  “I77F- 


(22) 


where  dx  * dx  i (Again,  dx  is  a negative  increment  in  the  integration) 

*3  » (XI  - x)  i + (Y1  + sR)  j + Zl  k 
and  | r^|  - (XI  - x)2  + (Y1  + sR)2  + Zl2 

Evaluation  of  the  cross-product  gives  dx  X r3  » Zl  dx  j + (Y1  + sR)  dx  k. 
Substituting  into  the  Integra’  Dives 


+ +rn  ““  Zl  dx  j - (Yl+sR)dx  k 

q3(P)  - ~ / 


(23) 


Air  o [(Xl-x)2+(Yl+sR)2+Zl2]3/2 

Evaluation  of  the  integral  is  identical  to  that  for  q^(P)  and  the  result 
in  component  form  is: 
y- component: 


-r0  Zl 


q3  (P)  j - — , . 

y 4x[(Yl+sR)2+Zl2] 


1 - 


XI 


z-component : 

q3  (P)  k 


+ro(Y3+sR) 

Air  [ (Yl+sR)  2+Zl2  ] 


[X12+(Y1+sR)2+Z12]1/2 


XI 

[X12+(Y1+sR)2+Z12]1/2 


j 


(24) 


(25) 


Total  Induced  Velocity 

In  summary,  the  total  induced  velocity  at  the  point  P is  the  vector 
sum  of  the  induced  velocities  created  by  the  bound  vortex  filament  and 

-V  ^ 

each  of  the  two  trailing  vortex  filaments,  that  is  q (P)  ■ q^(P)  + qj(P)  + 
(P) . The  total  velocity  at  point  P is  the  vector  sum  of  the  induced 
velocity  q(P)  and  the  free  stream  velocity  U*  ■ u.i  + v^j  + vjt.  If  the 
free  stream  velocity  is  directed  along  the  negative  x-axis  then  ■ - U^i 

■f 

and  the  relative  velocity  at  the  point  P is  UB  + q (P) . In  component  form 
the  relative  velocity  at  the  point  P is: 
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x-component:  u(P)  i « - U.i  + qj  (P)  i 

x 


rnp,  ,s  (s2  - y2)2  dy 
u(P)  - - U 4 -2-i  / 1 

Airs  -s  [Xl2+(Yl-y)24-Zl2]3/2 


y-component:  v(P)  J - q~  (P)  J + q,  (P)  j 

y Jy 


4-r0  zi  xi 

v(P)  1 - 

4n [ (Yl-sR) 2+Zl2 J [X12+(Y1-sR)2+Z12 ] 1/2 


r0  Zi 


” 1 

4ir  ( (Y14-SR)  24-Zl2  ] (X124-(Y1+sR)  24-Zl2  ] 1/;! 

*-component:  w(P)  k ■ q^(P)  k + q2  (P)  k 4-  q3  (P)  k 


r0  XI  s (s2  - y2)2  dy 
4 irs  -s  [Xl24(Yl-y)24-Zl2]3/2 


r0(Yl-sR) 

4tt[(Y1-sR)2+Z12] 


r0(Yl+sR) 

4ir[  (Y14sR)24Z12] 


(X124(Y1-sR)2+Z12]1/2 


(X12+(Y14-sR)24Z12]1/2 


Analysis  of  Induced  Velocity 

A computer  program  was  written  to  calculate  the  relative  velocity 
of  point  P as  the  point  was  moved  relative  to  the  horseshoe  vortex. 

Appendix  A contains  the  computer  program  which  was  used  to  calculate  the 
induced  velocity  from  the  equations  presented  in  the  preceeding  section. 

The  induced  velocity  q(P)  generated  by  the  horseshoe  vortex  model  is  plotted 

in  Fig.  9,  p.  25  » for  ■ 300.0  knots  and  in  Fig.  10,  p 26,  for  Uw  - 
200.0  knots. 

The  variation  of  induced  velocity  plots  were  generated  by  fixing  a 


value  of  XI,  ZI  and  then  computing  the  induced  velocity  for  the  range  of 
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P/SPAN  IS  BEHIND  AND  BELOW  HOTHERSHIP  WING. 

* l P/SPAN)  RATIO  = INF 

0 (P/SPAN J RATIO  = 1.61086 

0 l P/SPflN)  RATIO  = 1.09611 

A ( P/SPflN ) RATIO  = 0.88532 

X l P/SPAN ) RATIO  = 0.67452 

$ (P/SPAN)  RATIO  = 0.38822 

X l P/SPAN ) RATIO  = 0.219871  DOCKING  POINT) 


Figure  9.  Induced  Velocity  Distribution,  U ■ 300  Knots 
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values  of  Y1  shown.  The  Induced  velocity  was  non-dimensionalized  by  the 
free  stream  velocity  and  the  distance  Yl  from  the  plane  of  symmetry  of  the 
mothership  was  non-dimensionalized  with  the  mothership  wing  span. 

The  plots  o i non-dimensional  induced  velocity  verse  non-dimensional 
Y1  show  a variation  in  curvature.  The  variation  in  curvature  of  the 
velocity  was  caused  by  the  variation  in  the  effect  of  the  u,  v,  and  w 
components  of  induced  velocity.  The  u-component  comes  only  from  the 
elliptic  distribution  of  circulation  on  the  bound  vortex.  The  v-component 
comes  only  from  the  trailing  vortices  but  as  XI  increases  in  negative 
value,  the  fraction  XI  / [XI2  + (Y1  + sR)2  + Zl2]^2  increases  in  negative 
value.  The  fraction  is  subtracted  from  1 thus  reducing  the  v-component 
for  fixed  values  of  Yl.  Conversely  as  the  negative  value  of  XI  approaches 
zero  the  fraction  becomes  a progressively  smaller  negative  number  which 
is  subtracted  from  1 thus  increasing  the  v-component.  For  fixed  values 
of  XI,  Z1  the  v-component  of  induced  velocity  increases  as  the  trailing 
vortex  is  approached.  The  w-component  of  induced  velocity  comes  from 
both  the  bound  vortex  and  the  trailing  vortices.  The  merging  of  the  plots 
for  P/span  - 0.38822  and  0.21987,  at  values  of  Yl/mothership  wiug  span 
greater  than  + 0.20  is  the  result  of  the  v and  w components  magnitude 
dominating  the  u components  magnitude. 

At  values  of  P/span  - 0.67452  and  greater  the  relatively  flat  velocity 
variation  for  increasing  values  of  Yl  supports  the  assumption  of  constant 
velocity  across  the  RPV  span.  But  as  shown  in  both  velocity  plots  this 
assumption  breaks  down  as  the  docking  point  is  approached.  Variation  of 
velocity  across  the  RPV  span  becomes  important  if  the  RPV  is  out  of  the 
plane  of  symmetry  close  to  the  docking  point. 

The  induced  velocity  generated  at  ■ -200.0  knots  is  greater  than 
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IV.  Steady  State' Analysis  of  Aerodynamic  Interference 

The  flow  field  seen  by  the  RPV  during  approach  to  the  docking  point 
was  modeled  by  replacing  the  mothership  with  a horseshoe  vortex  as  pre- 
viously described.  It  was  shown  in  the  section  Analysis  of  Induced 
Velocity,  p.  24,  that  the  flow  field  is  essentially  constant  across  the 
RPV  wing  span  when  the  RPV  flight  path  is  within  a distance  of  + 0.1  of 
the  mothership  wing  span  from  the  mothership  plane  of  symmetry.  In 
analyzing  the  aerodynamic  interference  of  the  mothership  on  the  RPV  it 
was  assumed  that  the  RPV  flight  path  remained  in  the  plane  of  symmetry 
of  the  mothership. 

In  the  analysis  it  was  assumed  that  both  the  mothership  and  the  RPV 
were  flying  steadily  through  a stationary  atmosphere  (no  wind).  A trans- 
formation was  made  to  bring  the  aircraft  to  rest  and  give  the  air  a nega- 
tive velocity  equal  in  magnitude  and  opposite  in  direction  to  the  flight 
velocity. 

The  mothership  effects  a change  in  the  velocity  vector  seen  by  the 
RPV.  Both  the  magnitude  and  direction  of  the  velocity  vector  are  changed. 
To  show  frhe  effect  of  the  changed  velocity  vector  the  RPV  lift,  lift  in- 
duced drag,  and  aerodynamic  pitching  moment  were  computed  as  the  RPV 
approached,  the  docking  point.  In  the  analysis,  all  values  were  calcu- 
lated using  the  RPV  parameters  at  the  plane  of  symmetry  as  representative. 

The  approach  to  the  docking  point  that  was  found  to  be  best  during 
the  flight  tests  of  the  XF-85  is  shown  in  Fig.  11,  p.  30  (Ref  9:46-49). 
The  first  segment  of  the  flight  path  was  level,  about  5 times  the  diameter 
of  the  mothership  below  the  docking  point.  A gradual  increase  in  RPV 
altitude  was  then  made  to  bring  the  RPV  to  a point  approximately  one 
mothership  wing  mean  aerodynamic  chord  behind  and  level  with  the  docking 
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Figure  11.  RPV  Flight  Path 
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point.  The  RPV  was  then  flown  level  into  the  docking  point  at  a low 
closing  velocity  (less  than  5 knots).  For  the  analysis  the  RPV  approached 
the  dockiig  point  along  a level  flight  path  (0  ■ 0®  in  Fig.  11).  A level 
flight  path  was  chosen  to  remove  the  effects  of  climb  which  produced  dis- 
continuities in  the  data  at  the  points  where  the  flight  path  changed 
direction. 

A hypothetical  RPV  was  used  in  analyzing  the  variation  of  the  para- 
meters lift,  lift  induced  drag,  and  aerodynamic  pitching  moment  at  dis- 
crete points  along  the  flight  path.  At  the  discrete  points,  the  steady 
state  value  of  the  parameters  was  calculated.  Figure  12,  p.  32,  shows 
the  steady  state  equilibrium  forces  and  the  relevant  angles  for  the  RPV 
outside  and  inside  the  influence  of  the  mothership.  The  lift  of  the  RPV 
will  be  discussed  first,  followed  by  the  lift  induced  drag  and  aerodynamic 
pitch  moment. 

Variation  of  Lift 

The  lift  of  the  RPV  wing  is  given  by  L ■ l/2pV^  S C^.  The  coeffi- 
cient of  lift  is  a function  of  the  lift  curve  slope  (a)  and  the  angle 
of  attack  (a),  i.e.,  CL  ■ a(a-a^0).  For  the  special  case  studied  a 
symmetrical  airfoil  section  was  assumed  so  the  angle  of  attack  at  zero 
lift  (a^)  zero«  The  local  angle  of  attack  of  the  RPV  wing  changes 

from  the  free  stream  angle  of  attack  upon  approaching  the  mothership. 
Figure  13,  p.  33  defines  the  angles  and  their  relationship.  As  shown  in 
Fig.  13,  the  free  stream  angle  of  attack  (aB)  of  the  RPV  is  decreased  by 
an  amount  equal  to  the  angle  t Inside  the  influence  of  the  mothership; 
thus  the  local  angle  of  attack  in  the  mothership  influence  is  a ■ a.  - t. 

Variation  of  the  local  angle  of  attack  across  the  span  and  with 
position  relative  to  the  mothership  causes  a corresponding  change  in  the 
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a«  * Angle  of  Attack  of  the  RPV  in  the  Free  Stream 
Outside  the  Influence  of  the  Mothership. 

t ■ Angle  of  the  Relative  Velocity  to  the  Local 
Horizontal. 

o ■ Angle  of  Attack  of  the  RPV  with  the  Relative 

Velocity  Induced  by  the  Mothership,  (a  ■»  aa  - t) 


Figure  13.  Definition  of  Angles 
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coefficient  of  lift.  Since  the  lift  versus  a curve  slope  (CL  *a)  for 

a 

the  RPV  is  fixed  by  the  airfoil  section  used,  the  coefficient  of  lift 
varies  directly  with  the  change  in  local  angle  of  attack.  The  local  angle 
of  attack  has  been  assumed  constant  across  the  span  and  equal  to  the  local 
angle  of  attack  at  mid-span  so  for  the  ving 


C.  ■ aa  ■ a(a_  - t) 


(2J) 


The  lift  curve  slope  for  the  assumed  untwisted  elliptic  RPV  wing  is  com- 
puted using 


3Cj 

8a 


(30) 


1 + 


vAR 


where  a,,  is  the  theoretical  lift  curve  slope  and  AR  is  the  aspect  ratio 
(Ref  6:369).  The  value  of  a.,  is  2ir  (Ref  6:292). 

When  the  RPV  is  outside  the  influence  of  the  mothership  flying  in 
steady  level  flight  with  L ■ W the  lift  is  computed  from 


^ " 1/2  po  U»  S % “ 1/2  Po  U- 


Saa„ 


(31) 


Once  the  RPV  enters  the  influence  of  the  mothership  the  lift  is  computed 
from 

2 2 

(32) 


■V  * 1/2  % VLl  5 CL.  ' 1/2  4l  S 3 - '> 


REL  Ly  ’ "o  REL 

The  above  equation  shows  that  the  lift  in  the  influence  of  the  mothership 
is  directly  dependent  on  the  relative  velocity  vector  magnitude  and 
direction.  To  show  the  variation  of  the  local  angle  of  attack  caused  by 
both  the  change  in  the  relative  velocity  and  the  direction  of  the  relative 
velocity,  a was  computed  from 

W cos  T 


1/2  pc  S a V 


REL 


a » 


(33) 


GAW/AE/74-2 
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The  horizontal  tail  surface  also  adds  to  the  total  lift  of  the 
aircraft.  The  tail  contribution  to  aircraft  lift  is 


Lt  cose  - sine  <34 ) 

where  c is  the  dovnvash  angle  from  the  RPV  wing.  The  lift  of  the  tall  is 
computed  from 


L - 1/2  p„  V*  S C. 

0 REL  t L. 


(35) 


where  St  is  tail  area,  is  the  relative  velocity  of  the  wing.  The 
tail  coefficient  of  lift  is  computed  for  a symmetric  airfoil  from 


%.  ’ *t  “t  (36) 

where  at  ■ 3CLt/3at  and  is  the  tail  angle  of  attack.  The  slope  of  the 
tail  lift  coefficient  verse  curve  (at)  includes  the  variation  of 
velocity  between  the  tail  (V1)  and  the  wing  (Vr^*  The  tai*  an£le 
attack  is 

- a + ifc  - c (37) 

and  will  be  negative  to  produce  negative  lift  for  the  tail.  Thus  the 
lift  of  the  tail  is  given  by 

1 ■ 1/2  "o  VLl  St  at  <°  + lt  - £)  (36) 

where  it  is  the  tail  incident  angle.  The  drag  of  the  tail  is  computed 
as  shown  in  the  following  section. 

The  total  lift  of  the  RPV  is  the  sum  of  the  wing  generated  lift  and 
the  components  of  the  tail  forces  perpendicular  to  the  velocity  vector. 
The  total  lift  was  computed  by  the  formula 


L ,■  L + L cose  - D sine 
Total  W t t 


(3S) 
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Variation  of  Drag 

2 

Drag,  on  ti.e  RPV  wing  is  given  by  D ■ 1/2  p V S C , The  co- 

0 KhL  D 

efficient  of  drag  for  the  wing  is  given  by 

CD  " CD0  + k ^ WO) 

CpQ  is  the  coefficient  of  drag  at  zero  lift  and  k is  a factor  to  correct 
for  aspect  ratio,  thickness  effects,  pressure  drag;  all  drag  which  varies 
with  angle  of  attack  or  lift.  The  value  of  k is  computed  from 


k - -i—  (41) 

teAR 

where  e is  the  Oswald  airplane  efficiency  factor  and  AR  is  the  aspect 
ratio  of  the  wing  (Ref  10:7-12). 

For  a constant  CD0  and  k,  the  coefficient  of  drag  varies  as  the 
square  of  the  lift  coefficient.  Upon  substituting  for  C L the  coefficient 
of  drag  becomes 

% " So  + k a2°2  (42) 


Since  Cog,  k»  and  a are  constant  for  the  RPV  the  coefficient  of  drag  varies 

2 

directly  with  the  square  of  the  local  angle  of  attack  (a  ) and  the  square 
of  the  relative  velocity  in  the  mothership  influence.  The  drag  of  the 
wing  was  co&DUted  using 


2 1 ? ? 

D * 1/2  p V S(C  + a a ) 

W 0 reL  D„  TteAR 

0 


(43) 


The  drag  of  the  tail  is  computed  from 


D 


t 


1/2  po  Vrel  St 


where 


(44) 


(45) 
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* 

Ihe  value  of  k£  is  computed  as  for  the  wing. 

The  total  drag  of  the  RPV  is  the  sum  of  the  wing  generated  drag  and 
the  components  of  the  tail  forces  parallel  to  the  velocity  vector.  The 
drag  of  the  fuselage  is  not  included  in  the  model.  The  total  drag  was 
computed  by  the  formula 

°T0TAL  " DW  * Lt  8ine  * Dt  cose  (A6) 


Variation  of  Pitch  Stiffness 

For  an  aircraft  to  be  statically  pitch  stable  in  steady  flight  the 
slope  of  the  pitch  moment  coefficient  versus  angle  of  attack  curve  must 
be  negative  (C  < 0) . To  determine  the  pitching  moments  acting  on 
the  aircraft  are  determined.  Figure  14,  p.  38  , schematically  shows  the 
parameters  used  in  developing  the  pitch  moment  equation.  The  following 
development  uses  the  notation  found  in  Chapter  6 of  Ref  11. 

The  moments  are  taken  about  the  aircraft  c.g.  A positive  pitch 
moment  is  directed  parallel  to  the  positive  y-axis  and  causes  a nose  up 
pitch  angle.  The  body  contribution  to  the  pitch  moment  is  neglected  in 
this  analysis.  The  wing  contribution  to  the  pitch  moment  is 


Mjy  ■ Mac  + (Lw  cosa  + Dy  sino)  (h  - h^  )c 
w W w 

+ (Ly  sina  - Dw  cosa)z 

The  tail  contribution  to  the  pitch  moment  is 


(48) 


M ■ M - 1 [L  cos(a  - e)  + D sin(a  - e ) ] 
t ac„  t t t 

C (49) 

- cos(a  - c)  - Lt  sin(a  - c) 

The  moment  arm  lfc  for  the  tail  is  variable  because  the  c.g.  location 
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Figure  14.  Definition  of  Pitch  Moment  Parameters 
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moves.  The  mean  aerodynamic  center  location  is  fixed  for  fixed  wing 


aircraft  so  the  moment  arm  lt  can  be  written  in  terms  of  lfc  and  h as 


lt  - lt  - (h  - h )c  (50) 

w 

The  total  pitching  moment,  neglecting  body  contributions,  is 


M - \ + Mt 


(51) 


Thus, 


M ■ M + M + (L  cosa  + D sina](h  - h )c 
acw  act  1 w w ny' 


+ [L  sina  - D cosa]z 
w w 


- [L  cos  (a  - e)  + D sin(a  - e)]  l - (h  - h )c 
t t t nw 


- (D  cos(a  - c)  - L sin(a  - e)]z 
t t t 


The  pitching  moment  coefficient  (Cm)  is  obtained  from 


(52) 


H 


1/2  P vl  S c 


(53) 


To  obtain  the  slope  of  the  pitching  moment  (M)  versus  a curve  (C..  ) 

a 

the  pitch  moment  equation  is  differentiated  with  respect  to  a.  The  re- 
sult is 


3^aCy  ^act 

— . JL  + 


da  3a 


3a 


(54) 


3k,  3D 

+ { — — cosa  - L sina  + — — sina  + D„  cosaj(h  - h_  )c 
3a  w 3a  w nv 


3Lw  3Dw 

+ ( — - sina  + Ly  cosa - cosa  + D sinajz 

3a  3a 


39 


GAW/AE/74-2 


3Lt 

I — - cos  (a  - e)  - L 8in(a  - rWi  _ Jl  \ ^ 8Dt 
3a 


e)(l  ) + — - 8in(a  - e) 
da  3a 


+ Dr  cos  (a  - c)(l  -i£  )]{T  - (h  - h )c} 

3a  c n„ 


ou  t . 

* I 7 — cos  (a  - c)  - D sin  (a  - E)(l  - __  ) 

*J°  3a 


3L. 


3c 


sin  (a  - E)  - Lt  cos(a  - e)  (1 )]z 

3a  t 


The  partial  derivatives  — — — — — 


3D. 


follows: 


3a  ’3a  * IT  ’ IT  evaluated  as 


I . „2 


• \ - I p s CL  - ~ P uf  s 


aa 


3Lv  1 2 

r*rpli  Sa 
3a  2 ® 


(55) 


. n 1 2 
b*  Dw  " T P U SC 


D " 2 P U-  S(CD0  + k^a>2) 


5DW  2 2 

— 2 - P l£  k a a 

3a 


(56) 


c.  Lt  - i p U*  Se  ^ . i „ „2  st  ,t(.  + lt 


- t) 


3Lt  1 2 3c 

~--p  u„st  at (1  ) 

3a  2 c c 3a 


(57) 


d>  + a/, 

3D 


~ * P St  kt  At2(a  + i ~ c)(l  - — ) 

3a 


(58) 


The  total  equation  for  the  partial  derivative  of  the  aerodynamic  pitch 
moment  will  not  be  written  out  because  of  its  length.  The  value  of 


is  obtained  from 


Jma 
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C_  * 


1/2  p uj  S c 3a 
and  equations  54,  55  , 56  , 57,  and  58. 

In  the  computation  of  the  aerodynamic  pitch  moment  (C^)  and  pitch 
moment  verse  a curve  (Cj^)  the  free  stream  angle  of  attack  and  velocity 
were  used  to  obtain  the  free  stream  values.  To  obtain  the  values  of 
Cjj  and  in  the  influence  of  the  mothership  the  local  angle  of  attack 
and  relative  velocity  were  used.  A symmetrical  airfoil  section  was 
assumed  for  both  the  wing  and  tail.  The  non  uniform  flow  field  created 
by  the  mothership  has  the  effect  of  a cambering  the  symmetrical  airfoil. 
The  camber  created  by  the  non  uniform  flow  field  is  inversely  proportional 
to  the  radius  of  curvature  of  the  flow  field.  The  change  in  flow  field 
direction  was  small  which  implied  a large  radius  of  curvature.  For  this 
model  the  camber  effect  was  assumed  to  be  negligible  so  both  tne  moment 
about  the  aerodynamic  center  and  the  partial  derivative  of  moment  with 
respect  to  a were  assumed  to  be  zero. 


Analysis  of  Hypothetical  RPV 

To  study  the  change  J.n  aerodynamic  forceB  on  the  RPV  as  it  approached 
the  docking  point,  a hypothetical  RPV  was  assumed.  The  number  of  vehicle 
parameters  specified  for  the  analysis  was  kept  to  a minimum  and  equations 
were  used  in  the  calculation  of  the  additional  parameters.  The  parameters 
defining  the  RPV  used  in  this  analysis  are  given  in  Table  1,  p.  42  . The 
values  that  were  specified  are  set  apart  from  the  calculated  values.  The 
parameters  given  in  Table  1 were  assumed  without  performing  a preliminary 
design  of  the  RPV. 

A computer  program  was  written  to  calculate  the  aerodynamic  forces 


on  the  RPV.  The  free  stream  equilibrium  values  of  the  aerodynamic  forces 
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Table  1 


A.  MOTHERSHIP 


Input  Calculated 

o.  "5*  T0  - 861.0809  ft2/sec 

Aspect  Ratio  - 10.6 

Wing  Span  - 132.6  ft 

Weight  - 108,000  lb. 

U * 300  knots 

CO 

B.  RPV 


Input 

Wing  Span  ■ 20.0  ft 
Root  Chord,  Wing  ■ 4.0  ft 

Cjjq^  * 0.008 

ew  - °-8 
CMo  " °*° 

Root  Chord,  Tail  ■ 2.5  ft 

So, " °-008 

et  ■ 0.8 

1,  - 4.078* 

L 

VH  - 3.5 

Tt/c  - 3.5 

Z ■ 0.5  (Low  Wing) 

Zt  - 0.25 
h - 0.15 
h - 0.25 

QW 

K - 0.40 
t»WB 

Weight  - 6000  lb. 

Initial  Closing  Velocity  ■ 5 knots 


Calculated 


Wing 

Aspect 

Ratio  ■ 6.4 

Wing 

Area  - 

62.8  ft2 

Wing 

M.A.  Chord  » 3.4  ft 

Wing 

CL  * 
La 

a ■ 0.0835/deg 

Tail 

Aspect 

Ratio  ■ 4.2 

Tail 

Span  ■ 

8.2  ft 

Tail 

Area  * 

16.2  ft2 

Tail 

M.A.  Chord  - 2.12  ft 

Tail 

CL  " 

at  ■ 0.0743/deg 

Tt“ 

11.88 

ft 

r - 

0 

332.76 

ft^/sec 

«.  • 

3.87* 

Cl" 

0.324 

Cd- 

0.0146 

Cjj  - -0.320 


l! 

I !' 

; 


; 


I 
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were  calculated.  Then  the  equilibrium  values  for  the  aerodynamic  forces 
were  calculated  inside  the  mothership  influence.  To  show  the  variation 
in  values  along  the  assumed  flight  path,  the  point  of  analysis  was  moved 
in  discrete  steps  until  the  docking  point  was  reached.  The  last  point 
on  each  of  the  plots  in  Figures  15  through  24  correspond  to  the  docking 
point. 

In  the  analysis  the  RPV  was  moved  from  a point  approximately  1.3 
mothership  wing  spans  to  the  rear  and  level  with  the  docking  puint.  The 
gradual  increase  in  altitude  indicated  in  the  description  of  the  flight 
path  corresponds  to  a distance  equal  to  approximately  0.38  wing  spans. 

If  the  RPV  had  been  allowed  to  make  the  gradual  incraase  in  altitude  at 
a climb  angle  equal  to  the  deflection  of  the  velocity  vector,  a distance 
equal  to  approximately  43  wing  spans  would  have  been  required.  The  de- 
crease in  RPV  velocity  as  it  approaches  the  docking  point  is  shown  in 
Figure  15,  p.  44  . The  flight  tests  of  the  XF-85  showed  a decrease  in 
velocity  as  the  docking  point  was  approached.  Figures  16,  17,  and  18, 
p.  45  , 46  , and  47,  show  the  change  in  lift,  Cp,  and  angle  of  attack 
respectively.  The  increase  in  angle  of  attack  comes  primarily  from  the 
decrease  in  velocity  which  requires  that  the  angle  of  attack  be  increased 
to  maintain  the  necessary  lift.  The  increase  in  Cp  is  a direct  result  of 
the  increase  in  angle  of  attack.  The  increase  in  lift  is  the  combined 
result  of  the  decreased  velocity  and  increased  angle  of  attack. 

The  variation  of  Cp,  lift  induced  drag,  and  thrust  are  shown  in 
Figures  19,  20,  and  21,  p.  48— SQL  The  value  of  Cp  is  increased  as  the 
docking  point  is  approached  due  to  the  increase  in  angle  of  attack.  The 
lift  induced  drag  increases  as  the  docking  point  is  approached.  The 
increase  in  thrust  is  just  equal  to  the  increase  in  drag  in  steady  flight. 


43 


in 

in 

i- 

V— 

tn 

o 

o 

CO 

z 

z 

o 

in 

* 

X 

o 

o 

• 

■ 

• 

o 

300 

in 

II  II  II 

z 

a 

o_  >- 

in  I- 


q_  u 


89*101 

*-0I« 


09*101  SS  * 1 01  H*10l  98*101  82*10 

Ai I 0Q13A  d I HSd3H10W  / A1I3013A 


.00  0.40  0.80  1.20  1.60  2.00  2.40 

RPV  LOCRT I ON  / MOTHERSHIP  WING  SPRN 


.00  0.40  0.00  l.ZC  1.60  Z . 00  Z.40 

RPV  LOCATION  / MOTHERSHIP  WING  SPAN 


46 


0.40  0 . BO  l.ZO  1.60  Z.t 

LOCATION  / MOTHERSHIP  WING  SPAN 


RPV  SPAN  / MOTHERSHIP  SPAN  = 0. 1S083 

MOTHERSHIP  VELOCITY  = 300.0  KNOTS 

INITIAL  CLOSING  VELOCITY  = S.O  KNOTS 


RPV  LOCATION  / MOTHERSHIP  WING  SPAN 


P 


IIUILIJIII4J  !U.JIIUJ  il 

GAW/AE/74-2 


Ul  PUPW^IJIJ  I ' 


The  XF-85  flight  tests  showed  a drastic  increase  in  drag  and  thus  thrust 
as  the  docking  point  was  approached. 

The  variation  in  static  pitch  stability  (C^)  is  shown  in  Figure  22, 
p.  52  • As  the  docking  point  is  approached  the  slope  of  the  verse  a 
curve  decreases  in  negative  value  which  corresponds  to  a decrease  in 
aerodynamic  pitch  stability.  The  flight  tests  of  the  XF-85  indicated 
that  the  vehicle  became  marginally  stable  as  the  docking  point  was 
approached  due  to  the  excessive  control  manipulations  required  to  maintain 
a steady  position  for  the  final  hook-on. 

The  decrease  in  the  closing  velocity  is  shown  in  Figure  23,  p.  53  . 
The  closing  velocity  is  the  difference  of  the  x-component  of  the  RPV 
velocity  and  the  mothership  velocity.  The  time  required  to  reach  the 
docking  point  is  shown  in  Figure  24,  p.  54  . The  tine  plot  indicates 
the  time  span  covered  in  the  analysis. 


51 


RPV  SPAN  / MOTHERSHIP  SPAN  = 0.15083 

MOTHERSHIP  VELOCITY  = 300.0  KNOTS 


-ic  Pitch  Stability 


GAW/AE/74-2 


4 


I 


r 

f 

# 

i 

i 


i 


i 


i 


in 

to 

H* 

*— 

CO 

O 

o 

CO 

z 

z 

o 

X 

x: 

in 

•— « 

o 

o 

m 

m 

« 

o 

o 

in 

o 

CO 

II 

II 

ii 

z 

ac 

o_ 

>- 

to 

j— 

•— « 

Q_ 

o 

•— « 

o 

X 

>- 

_l 

CO 

I— 

UJ 

X 

1 — 1 

> 

UJ 

o 

X 

o 

t!> 

H 

—I 

z 

o 

UJ 

1 — 1 

X 

> 

to 

o 

X 

CL. 

_J 

M 

o 

z 

X 

OC 

in 

Q_ 

DC. 

X 

lO 

UJ 

1 — I 

X 

H- 

> 

•— I 

a. 

o 

z 

OL 

X 

1 — I 

00  * S2  01T02  00  * B l 00*01  oo:s 

f 03S  ) M30Q  01  3IJI1 


00*0° 


c 

i 


.00  0.40  0.80  l.ZO  1.60  2.00  2.40 

RPV  LOCATION  / MOTHERSHIP  WING  SPAN 


GAW/AE/74-2 


4 

V,  Conclusions  and  Recommendations 


Conclusions 

A steady  state  aerodynamic  analysis  was  performed  on  the  integrated 
launch/rccovery  system  proposed  for  in-flight  capture  of  a RPV.  The 
analysis  was  performed  by  modeling  the  mothership,  wing  only,  with  a 
bound  vortex  having  an  elliptic  distribution  of  circulation  across  the 
span  and  two  trailing  vortices  of  constant  circulation.  The  mothership 
flow  field,  which  was  the  induced  velocity  generated  by  the  mothership 
model,  was  used  to  calculate  the  steady  level  flight  aerodynamic  forces 
of  lift,  lift  induced  drag,  and  pitch  moment. 

The  comparison  of  results  obtained  from  the  analysis,  with  the 
XF-85  flight  tests  showed  that  the  trend  of  the  change  in  aerodynamic 
forces  was  correctly  predicted  for  the  lift,  drag,  and  pitch  moment. 

The  aerodynamic  analysis  was  performed  on  a hypothetical  RPV.  The  analysis 
may  be  applied  to  an  actual  RPV/mothership  system  to  obtain  an  initial 
understanding  of  the  variation  in  aerodynamic  forces  as  the  RPV  approaches 
the  mothership. 

Recommendations 

The  purpose  of  modeling  the  mothership  was  to  investigate  the  mother- 
ship  effects  on  the  RPV  as  it  approached  the  docking  point.  By  modeling 
the  mothership  as  a horseshoe  vortex  representing  the  wing,  the  only  effect 
studied  was  that  of  the  wing.  An  improvement  in  the  model  of  the  mother- 
ship  may  be  made  by  representing  the  body  with  a distribution  of  point 
sources  along  the  axis  of  the  fuselage  and  the  tail  with  a horseshoe  vortex. 
In  formulating  the  horseshoe  vortex  model  it  was  assumed  that  the  RPV  would 
be  farther  than  two  chord  lengths  away  from  the  center  of  pressure  of  the 
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wing.  This  may  not  be  the  case  for  some  potential  mothership  capture 
systems.  Thus  if  the  trapeze  docking  head  is  located  closer  than  two 
chord  lengths  from  the  center  of  pressure,  the  wing  probably  should  be 
modeled  as  a vortex  sheet.  The  horseshoe  vortex  model  of  the  mothership 
may  be  used  to  generate  a unsymmetrical  velocity  field  relative  to  the 
mothership  by  moving  the  point  of  analysis  off  the  plane  of  symmetry. 

The  unsymmetrical  flow  field  generated  may  then  be  used  in  analyzing  the 
dynamic  stability  of  the  RPV  as  it  approached  the  mothership. 
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Appendix  A 


Computer  Program 


A computer  program  was  written  to  perform  the  calculation  of  steady 
aerodynamic  forces  at  discrete  points  relative  to  the  mothership  model. 
Comment  cards  were  used  throughout  the  computer  program  to  indicate  what 
was  done  in  each  part  of  the  program.  A simplified  flow  diagram  is  given 
in  Figure  25.  A copy  of  the  computer  program  begins  on  page  60. 
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Simplified  Flow  Diagram  of 


Computer  Program  Model 


Input  Parameters 


Conversion  of 
Input  Parameters 


Compute  Mothership 
Flight  Parameters 


Compute  RPV  Vehicle 
Parameters 


Free  Stream  RPV  Level 
Flight  Calculations 


Compute  Velocity  Field 
Induced  by  Mothership 
as  RPV  Moves  Toward 
Docking  Point 


Compute  Variation  of  Steady 
Aerodynamic  Forces  on  the 
RPV  in  the  Influence  of 
Mothership.  Level  Flight. 


Plot  Aerodynamic  Forces 
Verse  Position 


Figure  25.  Simplified  Flow  Diagram 
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